Perturbation theory for sextic doubly anharmonic oscillator 
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A simple method for the calculation of higher orders of the logarithmic perturba- 
tion theory for bound states of the spherical anharmonic oscillator is developed. The 
structure of the perturbation series for energy eigenvalues of the sextic doubly anhar- 
monic oscillator is investigated. The recursion technique for deriving renormalized 
perturbation expansions is offered. 



I. INTRODUCTION 

The problem of the quantum anharmonic oscillator has been the subject of much dis- 
cussion for decades because of its importance in quant-field theory, molecular physics, and 
solid-state and statistical physics. However, while very extensive literature has been devoted 
to the one- dimensional consideration, relatively less attention has been given to physically 
more interesting but more complicated three-dimensional case. This applies to the doubly 
anharmonic oscillator of the type 

V(r) = -u 2 r 2 + Ar 4 + fir 6 , fi > (1) 

as well. In the one-dimensional space, there has been some interest in this theory both in 
its own right , and in the corresponding scalar field theory in 2+1 dimensions 
Besides double perturbation expansions in terms of the two coupling constants A and /i 
% [5|, this problem has been studied using the theory of continued fraction [6|, too. 

In the three-dimensional space, only exact solutions for the case when the parameters 

nnnnn 

satisfy certain relations have been obtained by a number of authors |7|, |8|, l9l . llOl . 1 1 If ] and first 
four eigenvalues have been calculated with the Hill determinant approach |l2l |. 

The main tool for studying energy eigenvalues and eigenfunctions of the bound-state 
problem within the framework of the Schrodinger equation is the logarithmic perturbation 
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theory (for references see 0, [3]). But the nonconvergence of obtained expansions compels 
us to resort to modern procedures of summation of divergent series, such as various ver- 
sions of the renormalization technique 5j]. However, to provide a reasonable accuracy these 
procedures need to know high orders of the perturbation series. Unfortunately, with the 
standard approach to the logarithmic perturbation theory we can easily obtain the higher 
order corrections only for the ground states, because the description of radially excited states 
becomes extremely cumbersome. 

Recently, a new procedure based on the specific quantization conditions has been pro- 
posed to get series of the logarithmic perturbation theory yia the ^-expansions within the 



framework of the one-dimensional Schrodinger equation 13| . Avoiding disadvantage of the 
standard approach, this straightforward semiclassical procedure results in the new handy 
recursion formulae with the same simple form both for ground and excited states. More- 
over, these formulae can be easily applied to any renormalization scheme of improving the 
convergence of expansions 5]. 

The object of this paper is to extend the above mentioned formalism to the bound-state 
problem for the spherical anharmonic oscillator with its subsequent applying to the doubly 
anharmonic oscillator. 



II. THE METHOD 

At the beginning we shall concentrate our attention on the bound-state problem for a non- 
relativistic particle moving in a central potential of anharmonic oscillator of a general form 
admitted bounded eigenfunctions and having in consequence a discrete energy spectrum. 
This potential has a single simple minimum at the origin and is given by the symmetric 
function 

V(r) = -mu 2 r 2 + ^ Vi r 2i+2 . (2) 

i>l 

Then after separating the angular part, the reduced radial part of the Schrodinger equation 
takes the form 

- ^U"{r) + + V(r)) U(r) = EU{r). (3) 

With applying the substitution, C(r) = hU'(r)/U(r), accepted in the logarithmic per- 



turbation theory, we go over from the Schrodinger equation (JS]) to the Riccati one 

W\r) + C 2 (r) = - ^ + 2mV{r) - 2mE. (4) 

We attempt to solve it in a semiclassical manner with series expansions in the Planck 
constant 

oo oo 

E = Y, E k h k , C{r) = C k (r)h k . (5) 

k=0 k=0 

In the complex plane, the number of zeros N of a wave function inside the closed contour 
is defined by its logarithmic derivative, C(r), through the relation 

— <tc(r)dr = — lc k {r)dr = hN. (6) 

27T1 / ^— ' / 

fc=0 

This quantization condition is exact and is widely used for deriving higher-order correc- 



tions to the WKB-approach Q, [3] and the 1/iV-expansions 
one important point to note. Because the radial and orbital quantum numbers, n and I, 
correspondingly, are specific quantum notions, the quantization condition (jHJ) must be sup- 
plemented with a rule of achieving a classical limit for these quantities. It is this rule that 
stipulates the kind of the semiclassical approximation. 

In particular, within the framework of the WKB-method the passage to the classical limit 
is implemented using the rule 

h — > 0, n — > oo, / — > oo, tin = const, M = const, (7) 

whereas the 1/iV-expansion requires the condition [lc]. \v\. fisj] 

h — > 0, n — const, I — > oo, tin — >• 0, hi = const. (8) 

The semiclassical treatment of the logarithmic perturbation theory proved to involve the 
alternative possibility: 

h — > 0, n = const, I = const, hn — > 0, hi — > 0. (9) 

Let us consider the latter rule from the physical point of view. Since hi — > as h — ► 0, 
the centrifugal term, h 2 l (I + 1) /r 2 , has the second order in h and disappears in the classical 
limit that corresponds to falling a particle into the center. This means that a particle 
drops into the bottom of a potential well as h —>■ and its classical energy becomes E = 
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181 ] . There is, however, 



mmV(r) = 0. It appears from this that the series expansion in the Planck constant for the 
energy eigenvalues must now read as E = YlT=i Eifi k - 

Upon inserting the ^-expansions for E and C(r) into the Riccati equation (jlj) and col- 
lecting coefficients of equal powers of H, we obtain the following hierarchy of equations: 

C 2 (r) = 2 mV(r), 
C' (r)+2C (r)C 1 (r) = -2mE 1 , 
C[(r) + 2 C (r)C 2 (r) + C 2 x {r) = - 2 mE 2 , (10) 



C' k -i{r) + Ci{r)C k „i{r) = -2 mE k , k > 2. 

i=0 

In the case of ground states, the recurrence system at hand coincides with one derived 
by means of the standard technique and can be solved straightforwardly For excited states, 
however, it is necessary to take into account the nodes of the wave function, that we intend 
to do by making use of the quantization condition 

It should be stressed that our approach is quite distinguished from the WKB method not 
only in the rule of achieving a classical limit but also in the choice of a contour of integration 
in complex plane. With a view to elucidating the last difference let us now sketch out the 
WKB treatment of the bound-state problem for the spherical anharmonic oscillator. In 
the complex plane, because the potential is described by the symmetric function ([2]), this 
problem has two pairs of turning points, i.e. zeros of the classical momentum. Therefore 
we have two cuts between these points: in the region r > as well as in the region r < 0. 
In spite of only one cut lies in the physical region r > 0, the contour of integration in the 
WKB quantization condition has to encircle both cuts for the correct result for harmonic 

n 

oscillator to be obtained |19| . 

In our approach, when a particle is dropping into the bottom of a potential well, these 
four turning points are drawing nearer and, at last, are joining all together at the origin. 
Hence, all nodes of the wave function are now removed from both positive and negative sides 
of the real axis into the origin and our contour of integration must enclose only this point 
and no other singularities. Further, let us count the multiplicity of zero formed at r = 0. 
For the regular solution of the equation ([3]), the behaviour r l+l as r — > brings the value 
I + 1. The number of nodes of eigenfunction in the region r > is equal to radial quantum 



number n. But because the potential is a symmetric function the same number of zeros 
must be in the region r < 0, too. Then the total number of zeros inside the contour becomes 
equal to N = 2n + l + 1. 

Taking into account the first order in % of the right-hand side, the quantization condition 
([6]) is now rewritten as: 

1 ' Ci(r)dr = 2n + Z + 1, — <bC k (r)dr = 0, k^l. (11) 



2ti-i J ^ ' 2tti„ 
A subsequent application of the theorem of residues to the explicit form of functions 
Cfc(r) easily solves the problem of description of the radially excited states. With that end 
in view let us consider the system ( jTUl) and investigate the behaviour of the function C k (r). 
From the first equation it is seen that Co(r) can be written in the form 

/ 2 oo \V 2 

C (r) = -\2mV(r)] 1/2 = -mtur 1 + -V^r 2 M = rV C, ( r 2i , (12) 

\ i=\ / i=0 

where the minus sign is chosen from the boundary conditions and coefficients Cf are defined 
by parameters of the potential through the relations: 

C° = -mu, C° = ^— C° P CI P - 2mvj , i > 1. (13) 

From (1121) we recognize that Co(0) = and, consequently, the function C\(r) has a simple 
pole at the origin, while Cfc(r) has a pole of the order of (2k — 1). Hence it appears that 
Ck(r) can be represented by the Laurent series: 

oo 

C k (r) = r x - 2k J2 C-r 2 \ k>l. (14) 
i=o 

With definition of residues, this expansion permits us to express the quantization condition 
(fTTj) explicitly in terms of the coefficients C\ as 

C k k _ 1 = (2n + l + l)6 1;k . (15) 



It is this quantization condition that makes possible the common consideration of the 
ground and excited states and permits us to derive the simple recurrent formulae. 



The substitution of the series (1T31) and <HM into the system (TTDTl in the case i ^ k — 1 
yields the recursion relation for obtaining the Laurent coefficients 

fe-l i 



c 



(3-2A; + 2z)Cf- 1 + ^^^a, 
i=i p=o 



k-j 
i—p 



+2^C° p Cl p + l(l + l)5 2 , k 5 , 
P =i 



(16) 



If i = k — 1, by equating the explicit expression for C|Li to the quantization condition (fT5|) 
we arrive at the recursion formulae for the energy eigenvalues 

fc-i 

2mE k = -Cti "EE C v°tL P ■ (17) 
i=o p=0 

Thus, the problem of obtaining the energy eigenvalues and eigenfunctions for the bound- 
state problem for anharmonic oscillator can be considered solved. The equations (Tl6l) and 
(|17p have the same simple form both for ground and excited states and define a useful 
procedure of the successive calculation of higher orders of the logarithmic perturbation 
theory. 



III. PERTURBATION CORRECTIONS TO ENERGY EIGENVALUES 

On applying the recursion relations obtained, analytical expressions for first corrections 
to the energy eigenvalues of the sextic doubly anharmonic oscillator (TjQ) are found to be 



equal to: 



1 + 2N 



(3-2L + 6r/) A 



2 



4 a; 
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£3 = 1 * ^ ((-21 + 9L- 17r7) A 2 + (15-6L + 10r/) u//i) , 
£ 4 = ^ (333 + 11 L 2 - 3 L (67 + 86 rj) + 3 77 (347 + 125 77)) A 3 

- 6 (60 + 3 (-13 + L) L+ 175 77 - 42 L 77 + 55 r/ 2 ) cj 2 A//J, (18) 



l + 2iV 
128 a; 11 



x 



(30885 + 909 L 2 — 21 L (613 + 330 rj) + rj (49927 + 10689 77)) A 4 
- 4 (11220 + 393Z 2 - 6L (1011 + 475t?) + 77 (16342 + 3129r/)) w 2 A 2 /j 
+ 2 (3495 + 138 L 2 + 4538 77 + 786 t/ 2 - 30 L (63 + 26 77)) u 4 ^ , 

where N = 2 n + I + 1, 77 = N (N + 1), Z, = 1(1 + 1), m = 1. 

As it is seen, that obtained expansion is the expansion in powers of the coupling constants. 
It is also evident that for the energy eigenvalues, when k — 1, we readily have the oscillator 
approximation 

Ex = [2n + l + - j w. (19) 



2, 

As a check of the obtained formulae we calculate the energy eigenvalue for the ground 
state of the anharmonic potential 

V(r) = ar 2 + h -r i + lr Q (20) 
3 9 



admitted the quasi-exact solution [11J. In this case the expression for first corrections to the 



energy eigenvalue with: n = 0, / = 1 take the form: 
E> 356 



2 " 24a ' 



E 



35 (-5b 2 + 6ac) 



3 



= ^ ^ (21) 



^5 = 
-^6 = 

Under the constraints 



96 V2a 5 
35 (475 b 3 - 864 a be) 
6912 a 4 

-35 (27565 6 4 - 67488 a 6 2 c + 17688 a 2 c 2 ) 

110592 V20 11 
35 (1451815 b 5 - 4482360 ab 3 c + 2489328 a 2 b c 2 ) 
2654208 a 7 



b = V2~c (2a + -v 7 ^) 172 (22) 
3 



we arrive at the analytical expression (with h — m — 1) 

E ,i = ^(2a+^v / 2^) 1 / 2 (23) 

derived with the supersymmetric consideration in 11]. 

Just as it was expected, the obtained series ([TBI for the sextic doubly anharmonic oscil- 
lator is divergent that compels us to resort to modern procedures of summation of divergent 
series. One of the most common among them are various versions of the renormalization 
procedures intended to reorganize a given series into another one with better convergence 
properties. 

It was found that the proposed technique is easily adapted to any renormalization scheme. 
Now we consider the case of the renormalization of a frequency. For this purpose it is enough 
to think of the harmonic oscillator frequency incoming in the potential as a function of 
Plank's constant variable, with subsequent its expansion in an fi-series. However, for later 
use it is more convenient to take 



u 2 



(24) 



k=0 



The relations for the coefficients C\ then reads 




1 i_1 

— (£C° p Cl p -2mv i ),i>l, 




1 



m 2 u\5i, k + (3 - 2k + 2i)C\ 



•k-l 



(25) 



k—l i i 



+ E E + 2 E + + , * > i, 



i=i p=o p=i 



but the recursion system for the energy eigenvalues does not change. 

The coefficients io\ are defined by the chosen version of renormalization. In practice, the 
one-parameter schemes are usually used. They are obtained with the truncation of the series 



where uq is a trial frequency and an order in h of the remainder is determined by the 
anharmonic term in a potential. In the case of our potential (JT]), taking into account the 
order in h of the dimensionless parameter of the expansion, we must put k — l. 

It has been recognized that for three-dimensional sextic doubly anharmonic oscillator 
the schemes of both the minimal difference, En{ujq) = 0, and the minimal sensitivity, 
(d/duo)E N (u> ) = 0,[5j appear to give approximately an equal accuracy in vast range of 
values of the anharmonic parameters. 

Typical results of the calculation are presented in the Table where the sequences of 
the partial sums of N renormalized corrections to the energy eigenvalues of the spherical 
anharmonic oscillator with the potential V(r) = \r 2 + Ar 4 + fir 6 is compared with the 
results, -Enum, obtained by the numerical integration for the values A = /i, (h = m = uj = 1). 
The free parameter ujq in the one-parameter scheme (1261) is defined order by order from the 
equation of minimal sensitivity, (d/duj^Ejy (a; ) = 0, with choosing such the root that obeys 
the condition of the flattest extremum \2l\ . 

From the Table it is seen that the one-parameter renormalization procedure gives the 
accuracy which is quite sufficient for practical purposes. 

In conclusion, a new useful technique for deriving results of the logarithmic perturbation 
theory has been developed. Based upon the ^-expansions and suitable quantization con- 
ditions, new handy recursion relations for solving the bound-state problem for a spherical 
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2 , 2*fe 
U Q + U k fl , 



(26) 



TABLE I: The sequences of the partial sums of N renormalized corrections. 
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anharmonic oscillator have been obtained. These relations can be applied to excited states 
exactly in the same manner as to the ground states providing, in principle, the calculation of 
the perturbation corrections of large orders in the analytic or numerical form. Besides this 
remarkable advantage over the standard approach to the logarithmic perturbation theory, 
our method does not imply knowledge of the exact solution for zero approximation, which 
is obtained automatically. And at last, the recursion formulae at hand are easily adapted to 
the use of any renormalization scheme for improving the convergence of obtained series. 

The proposed technique has been applied for investigation of the bound-state problem 
of the three-dimensional sextic doubly anharmonic oscillator. It has been observed that 
the perturbation series of the logarithmic perturbation theory are divergent but the use of 
the one-parameter renormalization procedure gives the accuracy which is quite sufficient for 
practical purposes. 
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